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The electromagnetic modes possibly unstable in strongly magnetized plasmas are identified. The 
regime where this instability might stand out compared to the incoherent electron-cyclotron radi- 
ation is explored. These modes are relevant to the inertial confinement fusion and the gamma ray 
burst. 
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The magnetic field is ubiquitous in plasmas. The 
plasma interaction with the magnetic field often deter- 
mines its dynamical behavior |l|-t5|. In particular, the 
magnetic field could convert the electron kinetic energy 
into the collective photons, as in the magnetron or the 
gyrotron [G, 7] and the free electron laser (FEL) (8l-[l0|. 
In the FEL, a relativistic electron beam gets periodically 
accelerated amplifying coherent electromagnetic (E&M) 
waves. The resulting lasers have various applications. 
One natural question would be whether the electron ther- 
mal gyro-motion in the magnetic field, not the specifically 
designed population-inverted plasma state, can excite an 
analogous process. 

In this letter, we examine the instability of collec- 
tive E&M waves arising from the electron thermal gyro- 
motion. The analysis of the electron motion in the pres- 
ence of a strong magnetic field leads to a theoretical 
framework similar to that of the Landau damping [Tl| , 
which reveals that there could exist numerous unstable 
E&M modes in relativistic plasmas. Our theory is similar 
to the well-known electron maser theory but our 

theory deals with the long time limit while the previous 
researches mainly focus on the short time limit, which 
will be discussed in detail at the end of Sec. (4). We 
identify the regime where the coherent radiation from 
this instability might be more intense than the inco- 
herent cyclotron radiation. Various implications of our 
study on the astrophysical and laboratory plasmas are 
discussed, including the short gamma ray burst the 
non-inductive current drive and the soft x-ray gen- 

eration in the inertial confinement fusion p^ . 

This paper is organized as follows. In Sec. (2), the 
theory of the instability is developed based on the Lan- 
dau damping theory for the non-relativistic plasmas. In 
Sec. (3), the fully relativistic theory is developed. This 
section is the major result of this paper. In Sec. (4), we 
identify the instability inherent in the relativistic plas- 
mas. In Sec. (5), we discuss the implications of our re- 
sults. 



In this section, we consider the non-relativistic plasma. 
While we will conclude that the theory developed in this 
section is inadequate, the intuition derived from this de- 
velopment is useful for the fully relativistic theory in the 
next section. We start with a non-relativistic electron 
under the magnetic field Bq = Bqz. The equation of 
motion is 



''~dt 



(Eo + ^ X Bo) , 



where rrie (e) is the electron mass (charge), Eq — 0, 
V is the electron velocity, and c is the speed of light. 
The zeroth-order solution is v^\t) — VpCos{uJcet + 4>o), 
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where 



ujce = eBo/mgC, Vp is the constant perpendicular ve- 
locity, and (po is the initial phase angle. Now, as- 
sume that this electron interacts with a linearly polar- 
ized E&M wave propagating in the positive z-direction: 
Ex{z,t) — EiCOs{kz — (jJt), Ey = Ez = 0, By(z,t) = 
Ei{ck/uj) cos{kz — u!t), and = = 0. The first order 
linearized equation is 
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Expanding the momentum equation of the second order 
in the z-direction and averaging it over c/jq, we obtain 
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where a = kvoz — uj + uJce and /3|| = v^z^ /c. Eq. ([T]) 
is exact to the second order in Ei. The first and the 
second term of the right hand side are from (v*^°^/c) x B 
and the third term is from (v^^^/c) x Bq in Eq. (P). The 
resonance condition, a — Q, leads to w — kv^z — uJce, and 
the resonance electron velocity is Vr = {uj — uJce) /k. Note 
that Vr is positive (negative) if > {uJ < ujce)- 
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FIG. 1: An electron emits a photon in the positive z-direction 
as it transitions from 1 to 2. See more details in the text. 



Eq. ([TJ shares similarity with the Landau damping 
analysis for a Langmuir wave [ill, 13 • From Eq. ([T]), 
we obtain the kinetic energy loss rate of electrons over 
the electron distribution in the limit of > 1 llal: 



de 



pe 

di. ~ 2 4^~fc2" 



ck 



2 dv 



ck. (2) 



where = v^/c^ = (v^ + vl)/c^, ^|| = vjc, uj"^^ = 
Annee'^ /nie is the plasma frequency, / is the electron dis- 
tribution function with the normalization of J fd^v — 1, 
and is the ensemble average with Vz — Vr or (A) = 
J AS(v — Vr)d^v. By employing the full relativistic equa- 
tion |12| and taking the appropriate classical limit, we 
also verify that Eq. ^ is correct in the energy conser- 
vation if the exchange of the electron parallel and per- 
pendicular energy with the E&M wave is fully taken into 
account. It should be noted that this equation of the 
energy exchange is well-known in the maser theory [l^ . 
The case when Vr < also seems contradictory as the 
E&M wave (the electron) gains (gains) the momentum 
(the parallel kinetice energy). This can be resolved by 
considering the perpendicular electron motion. The elec- 
tron kinetic energy in the perpendicular direction acts 
as an energy storage, by which the energy difference be- 
tween the electron kinetic energy in the z-direction and 
the energy gain of the E&M mode can be accounted for. 
We demonstrate this idea by a single particle simulation 
(Fig. [5]). As the resonant interaction between the electron 
and the wave progresses, the electron loses the momen- 
tum in the z-direction and, at the same time, gains the 
kinetic energy in the same direction. However, simulta- 
neously, the electron loses more energy in the perpendic- 
ular direction. The ratio of the perpendicular energy loss 
to the parallel energy gain is roughly c/vr, which is con- 
sistent with the momentum and the energy relation in 
the quantum mechanics; While the electron gains the ki- 
netic energy in the z-direction, it loses the energy in the 
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FIG. 2: The energy of an electron in the presence of the 
magnetic field and the E&M wave. The magnetic field is 1 
T so that cjce = 0.025 THz, and to = ck = 1/1.1 ojce- The 
intensity of the E&M field is 2 MW/cm^ sec. An electron has 
the initial perpendicular energy 80 keV and parallel energy 
4.05 keV. The x-axis is the time normalized by (cfc)~^ and the 
y-axis is the change in the parallel, perpendicular and total 
energy of the electron in the unit of keV. In this particular 
simulation, the electron parallel kinetic energy increases while 
its total (perpendicular) kinetic energy decreases. 



perpendicular direction by transitioning from a higher 
energy Landau level to a lower one. 

Assuming the energy density of the wave is given as 
e^j = QE'^/Att, we arrive at the wave instability growth 
rate (Landau growth rate) by equation the kinetic energy 
loss rate to the wave growth rate, Fe^ — {de/dt): 



ck. (3) 



In the case of a Langmuir wave, no instabilities exist if 
the electron distribution is peaked aX v — Q and mono- 
tonically decreases with v because the wave always sees 
the negative slope of the electron distribution at the res- 
onance. However, this is no longer the case for the E&M 
mode. If a; < Wce, Vr and k can be of the opposite signs 
and an amplification of the E&M wave (as well as the 
damping) can occur. An E&M mode propagating right- 
ward interacts resonantly with the electrons of certain 
negative velocity (Fig. [ij, and extract the momentum 
from the resonant electrons. The second term on the 
right-hand side of Eq. ([3]) is always a damping term, but 
the first term could be an amplification or a damping 
depending on the sign of df /dv. 

If C is positive, Eq. © predicts that the E&M 
wave becomes unstable in a Maxwellian plasma with 
Te > meC^/2, seemingly contradicting the Gardner's 
constraint [13, [3 which states that a plasma of an 
isotropic and monotonically decreasing distribution is 
stable. This apparent contradiction can be resolved with 
a proper estimation of the free energy in the plasma. If 
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an E&M wave is present in a plasma, the plasma ac- 
quires some free energy from the wave oscillation, and 
the E&M wave can extract this free energy up to the 
maximum intensity imposed by the Gardner's constraint. 
More specifically, let us assume that the electron distri- 
bution is Maxwellian for t < 0, and that an external 
E&M wave of is suddenly turned on at t = 0. For 
t > 0, the electrons acquire an additional kinetic energy 
of Se = CE^/iTT - = (C - 0.5(1 + iu}/ck)^))El/4TT, 
where e„ = (l/87r)i;2(l + (uj/ckf) is the E&M field en- 
ergy density. Let us assume that the E&M wave grows 
from Ex to Ef^ due to the instability. Then, the free 
energy drained into the E&M field energy is given as 
Se^ = 0.5(1 + {Lo/ckf){E% - ED/Att. We obtain the 
maximum wave intensity from the condition 5ew < Se, 
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III. 



LANDAU THEORY: FULLY RELATIVISTIC 
ELECTRONS 



While Eq. ([3]) is mathematically correct, the instabil- 
ity condition for the Maxwellian plasma predicted by the 
theory is Tg > meC^/2, so that the relativistic consider- 
ation is necessary. The momentum equation for a rela- 
tivistic electron is 



d7v 

nip — ; — = — e 
dt 



E 



X B 



- e- X Bo, 
c 



(5) 



where = l — v^/c^ is the relativistic factor. Following 
the same steps as in the classical case for a linearly polar- 
ized E&M wave but retaining only the resonance term, 
the electron energy loss rate after the average over the 
initial phase is 
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where f7(cj, fc, /3||) = {c^k'^/uj){l - uj(3\\/ck) - ujceh{.P), 
and a = kvz — uj + LOceH- Eq. ^ is exact to the second 
order in E\. In the regime where the gyro- frequency is 
very high so that ai 3> 1, the following approximation 
can be used nm: 
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where is a smooth function of the velocity /3 = v/c, 
dg/da — (V^a • \7 pg/\\i'a{/3)\'^) is the derivative of the 
g in the direction of the gradient of a. Employing the 
energy conservation, the growth rate of the E&M wave 
can be obtained by equating the energy growth rate of 
the E&M wave to the kinetic energy dissipation rate [iB] 



\-^"leC )AVG, 
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where Qavg is the ensemble average over the electron 
distribution. For ai 3> 1, the growth rate F is given 
from Eqs. dfi]), and dH) as 16] 
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where /3 = v/c, and 5(/3) = — uj/ck + ujce/{ckj{f3)), 
f is the electron distribution with the normalization of 
/ fd^P = 1, and {A)s=o = jS{S)Ad^/3 is the integration 
of the velocity space with the constraint S' = 0. Depend- 
ing on whether more electrons lose or gain the energy in 
the resonance boundary S{f3) = 0, the growth rate, F 
in Eq.®, could be of either sign. Let us consider the 
classical limit of Eq. (|9]), S = /3|| — uj/ck — uJce/ck and 
7 = 1. The resonance boundary is the same as the case 
in Eq. and Eq. (0) is simplified to 
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where O = (c^fc^/a;)(l — a;/3||/cfc) — ujce- Eqs. © and 
(ITU|) are the major result of our analysis. While Eq. © 
is similar with Eq. ([3]), there is a crucial difference; at 
UJ = ck, there could be an instability from Eq. 1^, which 
is not the case with Eq. One of the key limitations 
of Eq. ^ is the fact that the change of the cyclotron 
resonance frequency due to the change in the relativistic 
electron mass was not properly taken into account. 

As the distribution function described by the Vlasov 
equation is incompressible in the canonical coordinate 
of {p,q), the Gardner's constraint is partially relevant 
for the relativistic Maxwellian plasmas. Consider an ini- 
tially isotropic Maxwellian plasma, which has no collec- 
tive waves at its initial state. In the absence of the E&M 
waves, the electron kinetic energy is a function of the 
canonical momentum, and, due to the Gardner's con- 
straint, an certain amplified E&M wave cannot escape 
from the plasma if the final plasma has no collective E&M 
wave. However, a steady-state plasma without any col- 
lective wave rarely exists. Let us consider a plasma ini- 
tially with some collective waves of appreciable energy 



(more specifically plasmons and photons with w = Wpe). 
Noting that the Vlasov equation is compressible in the co- 
ordinate of (r, v) and the electron kinetic energy is given 
as E — ^Jm^d^ + (?{p— eA(r, i)/c)2, the Gardner's con- 
straint no longer applies because the electron kinetic en- 
ergy is now a function of the wave vector. Some col- 
lective E&M waves could be radiated while some other 
waves still remain in the plasma. Eqs. ([9|) and pO| sug- 
gest that the collective E&M wave is a channel through 
which the plasma free energy can be drained quickly. 



IV. 



INSTABILITY 



One notable consequence of Eqs. © and (|10l) is that 
an unstable E&M mode might exist even in Maxwellian 
plasmas. First, consider the case cfc < a; in the classical 
limit given in Eq. (jlOp . assuming the electron distribu- 
tion is an isotropic Maxwellian with the electron thermal 
velocity vth. In contrast to the example given in Fig. [TJ 
the first term on the right-hand side of Eq. (ITUl) is pos- 



itive \i Vr — Uj/k — UJce/k 

growth rate is given as 



> 0. For a positive Vr, the 
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where Pth = I'th/c, /?,- = Vr/c, and /||(t;) = 
l/A/27r?;tft,exp(— w2/(2w^^)). The condition for the insta- 
bility (r > 0) is 



Q\/ck + ck/u) 



(12) 



For an instability to be appreciable, ,f{vr) should be 
appreciably large or Eq. ((T^ should be satisfied near 
13 J. = j3th ■ The maximum possible growth rate at Vj. — Vth 
is given roughly as 



0.19 X 
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where it is assumed that C, 1. For an anisotropic 
plasma, the instability condition is 



1 - ( _ckVr^ 
C J Ld C 



(14) 



where ?;|| (wj_) is the electron thermal velocity in the par- 
allel (perpendicular) direction. If v\\/v± <C 1, the insta- 
bility range becomes much wider than the isotropic case, 
and the maximum growth rate is given as 
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For the fully relativistic relation (Eq. ©), the stability 
analysis gets more complicated. Here we consider only 
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FIG. 3: The contours of the resonance boundary (S'(/3) = 0) of 
the fully relativistic plasma for the case of cj > cfc > cjce- (a) 
u)/{ck) — 1.1 and ujce/ick) decreases from 0.9 (the outer-most 
ellipse) to 0.5 (the inner-most one) by 0.1, and (b) ijjce/{ck) = 
0.9 and ui/{ck) increases from 1.04 (the outer-most ellipse) 
to 1.32 (the inner-most one) by 0.04. Dotted lines are unit 
circles. The abscissa is /3x and the ordinate is /^y . 



one case of w > cfc > cjce- The Maxwellian electron 
distribution is given as /(/3) = 7*^ exp(— 7/ T,), which is 
the so-called Maxwell- Juttner's distribution [1^. When 
7o — Tf.jm^^c? > 1, it is peaked at 7 = 870 with the width 
of 6(3 = s/jq. Here, we assume that, at the resonance, 
the derivative of the distribution / is large dominating 
other terms. Then, Eq. ^ can be simplified to 
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2 c^k^ 
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Assuming /3 = 1, the condition F > can be recasted as 
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When uj > ck > ujce, the resonance boundary is of an 
elliptical shape (Fig. [3]). The resonance surface S has the 
maximum at Pmax = 1/ a/I -I- (wce/cfc)^. The necessary 
condition for the existence of the resonance is oj/ck < 

I3max + l^ce/^max, where -/max = (1 " Pmax)'^^'^- If 

electron distribution has the relatively high slope near the 
resonance region {■jmax — iT^/meC^), Eq. ([TT]) is possible 
since jV^S"! ^ 1 near the resonance. The growth rate 
without the gyro-damping term can be estimated, using 
V^/ = 7'(3-A7)//3 and/3 = 1, as 



(cfc)2 



5{S) 



7^^1(3- A)/ 



dS/dp 



d'/3, (18) 



where (dS/d^) = (/3 • VS")/!/?]. 

A similar analysis with Eqs. ([T|) and ^ is reported 
in the context of the electron cyclotron maser instabil- 
ity [13; however, the major focus is in the short time 
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growth rate {at < 1) when the magnetic field is not in- 
tense and/or the electron density is low. In the iner- 
tial confinement fusion or the astrophysical dense plas- 
mas, the situation of the magnetic field of giga-gauss and 
the electron density exceeding lO^^cm"'^ is common, and 
the Landau damping analysis valid for at ^ 1 should 
be sought instead of the conventional approach [l^ ; 
Eqs. (jni) and pU)) provide a proper estimation of the col- 
lective instabilities of dense strongly magnetized plasma. 
While the instability analysis in this paper focused on 
Maxwellian plasmas, Eqs. ^ and pUj) apply to not only 
Maxwellian plasmas but also general relativistic strongly 
magnetized plasmas such as jets, shock regions and accre- 
tion disk where the non-Maxwellian electron distribution 
is common. 

An E&M mode from the background noise could grow 
with the rate given in Eq. (|3]). This radiation could be 
more explosive than the incoherent cyclotron radiation. 
The ratio of the cyclotron radiation power per electron 
P to the electron kinetic energy is Td = {2P/inVp) = 
(4/3)(fce^/meC^)w, where the electrons are assumed to be 
non-relativistic. Using Eq. ([T3| . we arrive at Tmax/^d — 
(ne/k^){^/u!). If Tmax/^ci > 1, the radiation from un- 
stable E&M modes could be more explosive than the cy- 
clotron radiation. 



V. CONCLUSION 

To summarize, the coherent E&M instabilities of 
strongly magnetized relativistic electrons are analyzed in 
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